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; Detailed Current Edit History 
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WIMaeS Incos COMPLEX SINE AND COSINE 
/1-002/ ; File: MTHCSINCO.MAR 


| Pan eenaeUnANnRbeNntbenenEneNnsenendednesenbenstebsnnnsaenentonsesesonnceunes 


COPYRIGHT (c) 1978, 1980, 1982, 1984 B 

DIGITAL EQUIPHENT EGORPORATION, ” MAYNARD. MASSACHUSETTS. 

ALL RIGHTS R 

THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED 
ONLY IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE COPYRIGHT N NOTICE. THIS sof Tw ARE OR ANY OTHER 
COPIES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
eoRPORATiCn. " NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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FACILITY: MATH LIBRARY 


ABSTRACT: 
Return the S$ 


1 f a complex number 
Return the CO 


NE o 
SINE * » complex number 


VERSION: 0 
HISTORY: 


AUTHOR: 
Jonathan M. Taylor, 19-JUL-77: Version 0 


MODIFIED BY: 


1008) NOS NISTORY °; Detailed Current acme "BSEb=ASBS PEMi2E ERTARTE. SRESRrMCSiNCo.maR:1 29° 3) 


| -SBTTL HISTORY ; Detailed Current Edit History 
dit History for Version 0 of MTHSCSINCO 

0000 

0000 


-001 - vedere cig number and copyright notice. The last edit 
002 = Add "*_" e PSECT directive. JBS 21-DEC-78 
number S a hh 0 was 3. JBS 16-NOV-7 
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B11 
AND COSINE 16-SEP-1984 01: AX/VMS Macro V04-00 oof 
=$EE=1382 te 44:38 :26 CMTHRTL.SRCIMTHCSINCO.MAR; 1 ree 3) 


.SBTTL DECLARATIONS 
: INCLUDE FILES: 
: OERR.MAR 


; EXTERNAL SYMBOLS: 
-GLOBL MT 


: PSECT D DE aR ATIO 


CL NS: 

-PSECT  MTRSCODE PIC, SHR, LONG, EXE, NOWRT 

: EQUATED SYMBOLS: 
argadr = 4 : offset from AP of arg addr 


+ OWN STORAGE: 
$ NONE 
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M005 MTHSCSIN = COMPLEX SINE etsy Orde: 38 MTHRTL.SRCIMTHCSINCO.MAR; 1 ° (4) Tat 
i ~SBTTL MTHSCSIN = COMPLEX SINE 
hg 
3 ; FUNCTIONAL DESCRIPTION: 
38 : MTHSCSIN computes the SINE of a COMPLEX number (r, i) as 
3¢ 3 result = (SIN(r) * COSH(i), COS(r) * SINH(i)) 
39 3; CALLING SEQUENCE: 
6 $3 § Sine.wfc.v = MTHSCSIN(arg.rfc.r) 
Be8 100 ; INPUT PARAMETERS: 
00 101 ; The one input parameter is the address of a COMPLEX number (r, i), 
4 19 3 where r and i are both single-precision floating point values. 
000 104 : IMPLICIT INPUTS: 
$000 198 3 NONE 
0000 106; 
0000 107 ; OUTPUT PARAMETERS: 
0000 198 3 NONE 
000 109; 
000 110 ; IMPLICIT OUTPUTS: 
000 111; NONE 
0000 M6 3 
44 113 ; COMPLETION CODES: 
000 114; NONE 
0000 115; 
9000 116 ; SIDE EFFECTS: 
0000 117 ; Signals: Reserved Operand if r or _i are invalid (-0.0) 
0000 Bt 3 MTH$_SINSIGLOS if iri > 2eP]*2##31, 
0000 119; Floating Overflow if i > 88.028. 
464 120 ; 
000 121 ;-- 
B88 1 ; 
O60 1 
OOFC 000 124 -ENTRY MTHSCSIN, “M<R2,R3,R4,R5,R6,R7> 
00000025'EF 16 0008 125 JSB WORKER : RO =SIN(r) 
00 1 § 3; R1 = COS(r) 
0008 1 3 R2 = SINH(i) 
sit 128 3 R3 = COSH(i) : 
50 3) 4G 00 129 R3, RO ; RO = SIN(r) * COSH(i) 
51 52 $e 3 } 9 MULF R2, R1 3; R1 = COS(r) * SINH(i) 
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M005 MTHSCCOS = COMPLEX COSINE 6-SEP-1984 94539238 EMTHRTL. SREIMTHCSINCO.MAR: 1 . a) 
: ! ; -SBTTL MTHSCCOS = COMPLEX COSINE 
. Ji y= 
; ! § 3; FUNCTIONAL DESCRIPTION: 
9 : ! 8 : MTHSCCOS computes the COSINE of COMPLEX number (r, i) as follows: 
O0r 129 3 result = (COS(r) * COSH(i), -SIN(r) * SINH(-i)) 
i 146 3; CALLING SEQUENCE: 
OOF 143; Cosine.wfc.v = MTHSCCOS (arg.rfc.r) 
OF 144; 
OF 145; 
OQOOF 138 3; INPUT PARAMETERS: 
QOOF 147 ; The one input parameter is the address of a COMPLEX number (r, i), 
sit 138 3 where r and i are both single-precision floating point values. 
QOOOF 150 ; IMPLICIT INPUTS: 
QOOOF 151; NONE 
OOOF 13 3 
OOOF 153 ; OUTPUT PARAMETERS: 
QOOOF 154; NONE 
QOOOF 155; 
OOOF 138 : IMPLICIT OUTPUTS: 
QOOOF 157; NONE 
QOOOF 158; 
OOOF 159 ; COMPLETION CODES: 
QOOOF 160; NONE 
OOOF 161; 
OOOF 166 3; SIDE EFFECTS: 
OOOF 163 ; Signals: Reserved Operand if r or i are invalid (-0.0) 
OOOF 164 ; MTHS_SINSIGLOS if iri > 2*P]*2e31, 
OOOF 165 ; Floating Overflow if i > 88.028. 
OOOF 166; 
QOOF 167 ;-- 
OOOF 168 
BaOF 192 
OOFC O00F 170 eENTRY MTHSCCOS, “M<R2,R3,R4,R5,R6,R7> 
00000025'EF 16 0011 171 JSB WORKER : RO =SIN(r) 
0017 \¥¢ 3; R1 = COS(r) 
Bai) 17 3 Ré = SINH(i) 
017 174 3 RS = COSH(i) 
ee 44 0017 175 MULF R1, R3 3 R3 = COS(r) * COSH(i) 
28 28 6 Bata 126 MNEGF RO, RY ; RO = -SIN(r) 
51 ° 4 1 a MULF3 Re R2, R1 3 R1 = -SIN(r) * SINH(i) 
50 65 dO B35, 178 MOVL R3, RO ; RO = COS(r) * COSH(i) 
04 0024 179 ET 
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-SBTTL WORKER = do all the work 


a Serer handler 


" MTHSFLAG_JACKET 
MOVAB G*MTHSSJACKET_HND, (FP) 


; set handler address to jacket 


; handler 
MOVL ergadr (Ar). RO 3; RO => (r, i) 
MOVF 4(RO), RO 3; RO = 7 
JSB MTHSEXP_R4 3; RO = EXP(i) 
DIVF3 RO, #1.0, R1 3 R1 = EXP(-i) 
SUBF3 R1, RO, RS 3 RS = EXP(i) = EXP(-i) 
MULF #0.5, R5 3 RS = (EXPCI) = EXP(-i))/2 
ADDF3 R1,_RO, R6 3 R6 = EXPC(i) + EXP(-i) 
MULF #0.5, R6 : RO = (EXPCI) + EXP(-i))/2 
MOVF pergeee ee. RO ; RO=er 
JSB MTHSCOS_R4 ; RO = COS(r) 
MOVL RO, R7 3: R7 = COS(r) 
MOVF porgede eee RO ; RO=r 
JSB MTHSSIN_RS ; RO = SIN(r) 
MOVL R7, R1 ; R1 = COS(r) 
MOVL R5, R 3 Re = SINH(i) 
MOVL R6, R 3 RS = COSH(i) 
RSB 
-END 


set up error handler 
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MTHSCSINCOS COMPLEX SINE AND COSINE i oi St 4 9 91:14:16 ales Oe Macro Vv04-00 Page 7 
Symbol table 6-SEP-1984 11:21:26 CMTHRTL.SRC at MTHCSINCO.MAR; 1 (6) 
ARGADR = 00000004 
MIHSSJACKET_HND eeeeeree x (01 
MTHSCCOS QOOOOOOF RG 1 
MTHSCOS_R4 eeereree G 0 
MTHSCSIN 00000000 RG 1 
MTHSEXP_R4 eeteeree G 0 
MTHSSIN_RG eetereee G 0 
WORKER 00000025 R 01 
does wmoeanesnmecens + 
‘ ; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
. ABS. 00000000 ( 0.) 00 ¢ 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 0000006D ¢ 109.) O1¢ 1.) PIC USR CON’ REL Gch SHR EXE RD NOWRT NOVEC LONG 
feccuscesesesonessesesoses 
; Performance indicators ! 
Phase Page faults CPU Time Elapsed Time 
Initialization 33 00:00:00.09 00:00:01.38 
pounens processing 126 00:00:00.60 00:00:03.69 
84 Bp 288: 0.72 BB 120: Be = 18 
Synbol table sort 0 0:00:00.00 00:00:00.00 
s2 51 00:00:00.59 00:00:01.94 
ae table output 80:09:00 64 00:00:00.37 
Psect synopsis output 00:00:00.02 ee Se 
Cross-reference output 00:00:00.00 00:00:00.00 
Assembler run totals 300 00:00:02.06 00:00:09.56 


os, working rs Limit was 900 pages. 
2 bytes : pages) of virtua nore were used to buffer the intermediate code. 
There were 10 pages of symbol table space at yeceeee to hold 8 non-local and 6 local symbols. 
source lines were read in Pass 1, producing 14 object records in Pass 2. 
1 page of virtual memory was used to define 1 macro. 


poowceeceeoccccescececcccen} 


Macro Library name Macros defined 
“S255SDUAZB:(SYSLIBISTARLET.MLB;2 0 tS .. 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHCSINCO/OBJ=OBJ$:MTHCSINCO MSRC$:MTHJACKE T/UPDATE=(ENHS:MTHJACKET) +MS 
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